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Abstract 

o 

^ We show that strongly coupled field theories with holographic gravity duals at fi- 

nite charge density and low temperatures can undergo de Haas - van Alphen quantum 

^ oscillations as a function of an external magnetic field. Exhibiting this effect requires 

computation of the one loop contribution of charged bulk fermions to the free energy. 
The one loop calculation is performed using a formula expressing determinants in black 
hole backgrounds as sums over quasinormal modes. At zero temperature, the periodic 
nonanalyticitics in the magnetic susceptibility as a function of the inverse magnetic field 
depend on the low energy scaling behavior of fcrmionic operators in the field theory, and 
arc found to be softer than in weakly coupled theories. We also obtain numerical and 
WKB results for the quasinormal modes of charged bosons in dyonic black hole back- 
grounds, finding evidence for nontrivial periodic behavior as a function of the magnetic 
field. 
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1 Introduction 



1.1 Beyond universality: one loop effects in holography 

One objective of applications of the holographic gauge/gravity correspondence [1] to con- 
densed matter physics is the characterisation of exotic states of matter. Recent works have 
begun to uncover a rich structure in strongly coupled theories with holographic gravity 
duals at finite charge density. Initial studies [2 [3] focused on hydrodynamic aspects at 
higher temperatures, while many interesting ground states have emerged in later studies at 
low temperatures. When probed with charged scalar operators these theories can exhibit 
low temperature instabilities towards superconducting phases |ll|5l|6]. Equally interesting 
are the cases in which the finite density theory admits gapless charged scalar excitations 
but no superconducting instability [7J. When probed with charged fermionic operators, 
the response functions of the theory appear to indicate an underlying Fermi surface with 
non-Landau liquid excitations [51 |9l [101 E] • 

The recent discoveries listed above lead to a seemingly paradoxical situation. The pres- 
ence or absence of superconducting instabilities and Fermi surfaces is sensitive to the charge 
and mass of matter fields in the gravitational bulk spacetime [71 [11]. Equivalently, it is 
sensitive to the charge and scaling dimensions of low dimensional operators in the field 
theory. In contrast, the thermoelectric equilibrium and response properties of the theory 
are completely independent of these fields: they are universally determined by the Einstein- 
Maxwell sector of the bulk action. The fact that quantities such as the shear viscosity over 
the entropy density [12] or the electrical conductivity over the charge susceptibility |13] 
are identical for many distinct strongly interacting theories has been considered a robust 
prediction of sorts of applied holography. However, it seems unphysical that, for instance, 
the frequency dependent electrical conductivity should be independent of whether or not 
the theory has a Fermi surface. The latter property depends on the matter content of the 
bulk theory while the former does not. 

We take the viewpoint that the emergence of a universal thermoelectric sector from 
theories with radically distinct bosonic or fermionic response is an artifact of the classical 
gravity ('large in field theory) limit. In the absence of symmetry breaking condensates 
or relevant perturbations of the theory, the minimal gravitational dual at finite density 
and temperature is the charged AdS-Reissner-Nordstrom black hole (see e.g. |!l4!). In this 
background only the metric and Maxwell fields are nontrivial and all matter fields vanish. 
The fact that the Einstein-Maxwell sector does not source the matter fields is why the 
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equilibrium and thermoelectric linear response of the theory can be universal and blind to 
the matter content. Beyond the classical limit, however, it is clear that this 'decoupling' 
cannot continue to hold. All matter fields will run in loops and modify the gravitational 
propagator while all charged matter fields will modify the electronic properties of the theory. 
These effects are obviously small in the large limit, yet they may lead to qualitatively new 
physics. Furthermore, in any putative 'real world' application of holographic techniques, 
the desired value of N is unlikely to be large. 

This paper will initiate the study of bulk one loop physics in applied holography. These 
are a more general set of '1/A' corrections that go beyond those captured by including 
higher derivative terms in the gravitational action (see e.g. |15| [T6l [T7\ [18]). In particular, 
we will be interested in effects that are not captured by local terms in an effective action, 
as they involve loops of light fields. 

Consideration of quantum effects in the bulk potentially opens the Pandora's box of 
quantum gravity]^ Control of the ultraviolet properties of the quantum theory will ulti- 
mately require embedding computations into a consistent theory of quantum gravity. This 
may appear at odds with the 'phenomenological' approach to applied holography (e.g. |14] ) 
in which one restricts attention to a minimal set of fields needed to capture the physics of 
interest. At the one loop level these problems do not arise. Indeed the one loop physics 
of quantum gravity was fruitfully explored well before the availability of ultraviolet finite 
string theories [19]. The technical point that makes this possible is that functional deter- 
minants can be computed up to the renormalisation of a finite number of local couplings 
in the classical gravitational action. In particular, the nonlocal effects of interest to us are 
insensitive to the ultraviolet completion of the theory. 

1.2 Quantum oscillations as a probe of exotic states of matter 

In this paper we will consider the bulk one loop correction to the free energy due to charged 

matter. Our primary objective is to study the free energy as a function of an external 

magnetic field. Magnetic fields are fundamental probes of matter at low temperatures. 

The quantum Hall effect and closely related de Haas-van Alphen quantum oscillations are 

examples of phenomena in which Landau level physics reveals important information about 

the finite density system, such as the presence of a Fermi surface. 

^The 'box' of Pandora's box is apparently a mistranslation of the Greek word 'pithos' which refers to a 
large jar, often human-sized. As well as sickness and toil, the opening of the jar was also said to unleash 
Hope onto humankind. 
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In recent years, experimental studies of quantum oscillations have had a profound impact 
on our understanding of a variety of correlated electron systems. In the hole-doped cuprates, 
the observation |20| [21} |22| [23l IMl [25l |26] of quantum oscillations with a period indicative 
of 'small' Fermi surfaces has shown that the 'large' Fermi surface Fermi liquid state at 
large doping must be strongly modified in the underdoped regime. In the electron-doped 
cuprates, quantum oscillations with both small and large periods have been observed |27j . 
separated by a presumed quantum phase transition. In these contexts, it appears of interest 
to catalog the states of matter which can exhibit the quantum oscillations, apart from the 
familiar Fermi liquid suspects. Possible examples include fermionic matter coupled to gauge 
fields, or non-superfluid states of bosons such as vortex liquids or 'Bose metals'. 

It is not yet clear to what extent we can interpret the finite density matter of the 
present gravity duals in terms of the concepts mentioned in the previous paragraph, but 
it is our hope that a study of quantum oscillations will advance our understanding of such 
issues. The classical bulk (large A^) free energy is not manifestly written as a sum over 
Landau levels, as we shall see. The one loop correction to the free energy, in contrast, will 
naturally appear in this form. It follows that quantities such as the magnetic susceptibility 
can be expected to show novel qualitative features at the one loop level that are not visible 
classically. 

In section |2] we review the computation of the low temperature magnetic susceptibility 
for free fermions and bosons with a finite chemical potential. The case of fermions leads 
to de Haas-van Alphen oscillations (at the low temperatures we consider, these can also 
be thought of as quantum Hall transitions). We then compute the leading order in large 
N magnetic susceptibility in strongly coupled theories with gravitational duals in section 
[sj with no indication of quantum oscillations. We go on to consider the (bulk) one loop 
magnetic susceptibility in the strongly coupled theory. This is done using a, new to our 
knowledge, expression for determinants in black hole backgrounds written as a sum over the 
quasinormal modes of the black hole. The formula is derived in [28 and allows us to use 
the recent analytic results of [llj on fermionic quasinormal modes. Quantum oscillations 
are seen to appear from the one loop contribution of fermionic fields. We find that the 
periodic delta functions characterising the free fermion susceptibility at zero temperatures 
are replaced by power law divergences at strong coupling. In section |4] we numerically ex- 
plore the quasinormal modes of charged bosonic fields, discussing the possibility of periodic 
oscillations due to bosons also. 
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2 Free theories: fermions and bosons 



To introduce some of the techniques and concepts we will use later, we first exhibit the 
de Haas-van Alphen quantum oscillations in a more familiar setting. We will consider the 
cases of free charged bosons and free charged fermions in 2+1 dimensions. We work in 
Euclidean signature (t = —ir) and place the theory in a background chemical potential n 
and magnetic field B 

A = i^dr + Bxdy . (1) 

We are interested in computing the free energy at low temperature (T//x 0) as the 
magnetic field is varied at fixed chemical potential. 

It is convenient to treat the case of bosons and fermions simultaneously. For this purpose 
we can start with the Euclidean action for a complex scalar boson: 



and the following action for fermions 



d'^x 



W • ((9 - iqA) + 



These two actions give the free energy 
= Ttrlog [-V^ + 



^ tr log 



-V^ ±qB 



(2) 



(3) 



(4) 



where = dfj, — iqAf^. The only important difference between the bosonic and fermionic 
cases is that the bosons are periodic in the thermal time circle whereas the fermions are an- 
tiperiodic. The extra term appearing in the fermionic case is the magnetic Zeeman splitting 
of the spin degeneracy. This term will not qualitatively affect the quantum oscillations. 

The traces in Q can be computed as a sum over eigenvalues of the Laplace operator. 
The eigenvalues are given by 



V^^> + m^^> ± qB^ 



(5) 



where the it term should be added for fermions and is absent for bosons. We will retain 
this notation in the remainder of this section. The eigenvalue spectrum can be determined 
exactly by separation of variables in this equation. Let 



(6) 
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where A; G M and the thermal frequencies are 



LOn = 2i:nT (bosons) , (7) 

= 27r(n + \)T (fermions) , (8) 



OJ. 



for n G Z. The Xi{x) satisfy 

- X'l + q^B^x^X^ = K^Xt , (9) 

where we shifted the x variable so that x = x — k/qB. This equation for Xi is just the 
Schrodinger equation for a harmonic oscillator and therefore 

K,= \qB\{2l+l), (10) 

with£ G Z+U{0}. The eigenfunctions are Hermite polynomials X^(x) = e~l^-^l^'^/^i^^(y^|qi?|x) 
Putting the above together leads to 

A = + 2\qB\{l +\±\)- {itOn - q^J-f ■ (H) 

We see that the eigenvalue A will be independent of the momentum k. This is the degeneracy 
of the Landau levels. We can now check that the degeneracy is in fact 

(12) 

where A is the 2 dimensional area of the sample. To see this suppose that we had a finite 
sample of size Lx x Ly. The allowed values for the momentum would he k = 2Tmy/Ly for 
Uy G Z"^U{0}. The shift x — > x — k/qB we noted above is possible provided that k/qB < Lx- 
This places an upper bound on Uy leading to (12). 



Taking into account the degeneracy (12) of the Landau levels, one can perform the 



sum over eigenvalues to obtain the standard expressions for the free energy of bosons and 
fermions. For future comparison we will express the result in the following form 




(13) 



The upper sign is for bosons and the lower for fermions. For fermions one should addi- 
tionally let —>■ I X^^±, separating out the spin up and down contributions. A divergent 
temperature independent constant proportional to q/j, has been neglected. In the above 
result 



z.{i)=q^i±Jm^ + 2\qB\{e+^±^), (14) 



These values of are to be thought of as complex frequencies which give A = upon 



analytic continuation z = iiOn of (11 ). That is to say, they are solutions to the equations of 
motion, and zeroes of 

A(z, l)=rr?^ 2|g5|(^ + i ± i) _ _ g^)2 . (15) 

Expressing the free energy as a sum over complex frequencies that give zero modes of the 
differential operator is the key step that we shall generalise below at strong coupling. 



The sum over ^ in (13) will diverge. This is a temperature independent divergence, as 



can be seen by rewriting 

"log (1 T e-(^^-^^)/^) + log (1 T e-(^^+'?^)/^' 
+ (T independent terms) . (16) 



T j;iog(lTe-^*(^)/^) = T 



The finite temperature sums over i are now manifestly convergent. We introduced the 
energy of the £th Landau level 



e, = ^/m2 + 2|gi3|(^+i±i). (17) 

One can use a renormalisation method, such as zeta function regularisation, to control the 
zero temperature sums over Landau levels. At this point we should also comment on the 
zero magnetic field limit. The i? — > limit is to be taken keeping 

2\qB\l = fixed as B ^ . (18) 

In this limit the sum over the Landau levels becomes an integral over momenta 

^max. /-fcrnax. 

\qB\ V - / kdk, (19) 
« JO 



with /cinax. related to £max. via (18). The difference between bosons and fermions due to 



Zeeman splitting drops out in this large Landau level limit. 



From (13) and (16) we can see the de Haas- van Alphen magnetic oscillations in the case 



of fermions (the lower sign in these two equations). Take the T — > limit of (16) with 
fermionic signs and with fi and B fixed. Whether or not a given term contributes to the 
sum over i in this limit depends on whether —S£zizqfj, is positive or negative. If it is negative, 
then the exponential in (16) diverges and the term gives a finite contribution. However, if 
it is positive, then the exponential goes to zero, the argument of the logarithm goes to one, 
and hence the total term goes to zero. Therefore we have 

lim = — ^(9/^ — £i)(){QfJ' — e^) + • • • • (fermions, qfj, > 0) (20) 

T — ^0 ^TT 



Here 9{x) is the Heaviside step function and is equal to 1 if x > and zero otherwise. The 
dots denote analytic terms. As before for fermions ^£ = ^ X^^±. We see that the free energy 
changes nonanalytically whenever one of the Zi,{l) changes sign, say by tuning the magnetic 



field B. Note that this can only occur for one of the signs in (14), depending on the sign of 



q/i. In (20) we have assumed for concreteness that g/i > 0. Of course, these nonanalyticities 
will get smoothed out at any finite temperature. The jumps in the derivative clearly occur 
whenever a Landau level crosses the fermi energy. To see the oscillations themselves we 
should differentiate twice to obtain the zero temperature magnetic susceptibility 



X 



1\2 



E ' j ^i^f^ - + • • • , (21) 



where dots denote terms without delta functions. We can see that the susceptibility x shows 
a strong response with period 



A f i^l - ''^ - (22) 



where Ap = nkp is the cross sectional area of the Fermi surface, with kp = Ep — re? 



In (20) the zero temperature free energy is piecewise linear in the chemical potential. 
It we compute the charge density via p = d^/dp, then we find that the charge density is 
piecewise constant, with finite jumps at specific values of the magnetic field. These are the 
integer quantum Hall phases. 

The boson system is quite different. The system is only stable if £i > \qn\. If \qn\ 
becomes larger than eq then either the charged particles or antiparticles will condense, at 



any temperature. Using (16) the expression (13) is rewritten in the more familiar form 



\qB\AT 



2tt 

I 



E [log (l - e-^"^-"^)/^) + log (l - e-(^^+'?'^)/^)] + 9\p^^ . (bosons) 



(23) 

This last equation is recognised as the free energy of a gas of free charged particles and 
antiparticles. There are no jumps in the derivative, instead $7 diverges if becomes equal 
to one of the ££. 

Assuming that the mass is sufficiently large compared to the chemical potential so that 
the system is stable, the zero temperature free energy may be computed by, for instance, 
zeta function regularising the sum over Landau levels. One obtains 



^|g5|3/2 / 1 1 



m 



2 
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where the Hurwitz zeta function is defined by analytic continuation of 



Ch{s,x) = 



n=0 



{x + n) 



(25) 



The susceptibihty obtained by differentiating this expression twice is shown in figure [T] 
There are clearly no oscillations of the type obtained for fermions. The values of the 
dimensionless susceptibility appearing in the plot are seen to be small. Note that the 



chemical potential does not appear in ( 24 ) , so there is no charge density. The susceptibility 



in the plot is purely due to vacuum fiuctuations 
-0.002 n — . — . — . — . — , — . — . — . — . 



-0.004 
-0.006 
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Figure 1: The zero temperature magnetic susceptibility for bosons as a function of the 
magnetic field. The expression plotted has been made dimensionless by dividing by the 
sample area and multiplying by the boson mass m. 

In condensed matter applications, the theory Eq. ^ describes the superconductor- 
insulator transition of charged bosons at integer filling in a periodic potential; for this case 



Eq. (24) describes the diagmagnetic response of the insulating phase. 



3 Strongly coupled theories with gravitational duals 

3.1 The normal state geometry and large free energy 

In the previous section we reviewed the computation of magnetic susceptibility for free the- 
ories of bosons and fermions at finite chemical potential. We will now study the magnetic 
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susceptibility of certain strongly coupled field theories, again with a finite chemical poten- 
tial. Specifically, we study field theories which have large gravitational duals described 
'holographically' by Einstein-Maxwell theory in one dimension higher than the field theory 
(see e.g. [H] for a motivation of this dual description). We work with 2-1-1 dimensional field 
theories and hence 3-1-1 dimensional gravitational duals. 

Recall that our motivation is twofold. Firstly, we would like to see if any novel features 
arise in the magnetic response for theories that are stable against superconducting instabil- 
ities at finite chemical potential, despite having massless charged bosons [7]. Secondly, we 
would like to see if the putative Fermi surfaces identified in fermion spectral functions in 
[HI [9l [ini [TT] manifest themselves in the expected way as quantum oscillations. 

In the absence of superconducting instabilities, the state of the field theory is dually 
described by a solution to Einstein-Maxwell theory. We are interested in thermodynamic 
properties and so we shall work in the Euclidean theory. The Euclidean action is 



'2^(^+:&) + 4^^' 



(26) 



Here F = dA is the electromagnetic field strength. The Einstein equations of motion are 
while the Maxwell equation is 

V^F^" = . (28) 

The normal state at a finite temperature, chemical potential and magnetic field is de- 
scribed by the dyonic black hole metric (see e.g. 

ds^ = ^ (f{r)dT^ + ^ + dx'dx' ] , (29) 



^2 \^ J^j. 

with 

n„.^.(^.^A,l^)(^g,^Atl^(^i.j . ,30) 

together with the gauge potential 



A = in 



r 

1 



In these expressions we introduced the dimensionless quantity 



dr + Bxdy. (31) 



7^ = ^ , (32) 



which is a measure of the relative strengths of the gravitational and Maxwell forces. For a 
given theory, this ratio will be fixed. Some values arising in Freund-Rubin compactifications 
of M theory are described in [7] . 
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The field theory dual to this background has chemical potential ^, magnetic field B and 
a temperature given by the Hawking temperature of the black hole 

1 



(33) 



Note that whereas the chemical potential /x and temperature T have mass dimension one 
in field theory, the background magnetic field has mass dimension two. The free energy is 
given by evaluating the on shell classical action (see e.g. [29\ E]) 



1 + 



rln"^ 3rlB^ 



2K?'r^ \ 7' 7^ 
where A is the spatial area. From the free energy one computes the charge density 

1 d^n 2L2 



P 



and the magnetisation density 



m 



A dn 



1 dVLQ 
'A~dB 



^2 r^j"^ 



2L2 r+B 

K? 7^ 



(34) 



(35) 



(36) 



In these expressions, should be thought of as a function of fi,B and T via (33). 

Using the above results, the magnetic susceptibility x = —d^^o is easily computed from 
( 34 1 and ( 33 ) . The zero temperature result is plotted in figure [2j 



-2 
-3 

-5 
-6 
-7 



Figure 2: The zero temperature magnetic susceptibility to leading order at large A'^ as a 
function of the magnetic field. The expression plotted has been made dimensionless by 
dividing by the sample area and multiplying by the chemical potential /i. 
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Figure [2] is the leading order large N limit of the magnetic susceptibility]^ The plot 
is disturbingly similar to that for free bosons in figure [T| Note however that the strongly 
coupled theory is scale invariant, and so the only scale at zero temperature is the chemical 
potential fi, whereas in the free theory of the previous section we had a mass scale m. 

3.2 One loop (1/A^) corrections to the free energy 



The leading order at large result for the free energy, (34), clearly does not show any 
nonanalytic structure as a function of the magnetic field at low temperature. The magnetic 
susceptibility is correspondingly uneventful as shown in figure |2j 

We will show in the remainder of this paper that this uneventfulness is an artifact of 
the large N limit. As we mentioned in the introduction, a similar issue is known to arise in 
linear response. While the bulk Einstein-Maxwell theory captures all of the leading order in 
N electromagnetic and thermal response of the field theory, it appears to be independent of 
the spectrum (charges and scaling dimensions) of low lying fermionic and bosonic operators 
in the theory. Yet it is precisely this spectrum that determines whether or not there is 
a superconducting instability [7] and whether or not the fermionic response shows Fermi 
surface-like features A natural resolution to this tension is found in the fact that 

the Einstein- Maxwell and matter fields (fermions and bosons) are coupled in the bulk at a 
nonlinear level. Thus at higher orders in the expansion, or in higher point correlators, 
the matter fields will explicitly influence thermoelectric response. 

In what follows we consider corrections to equilibrium thermodynamic quantities, 
in particular the magnetic susceptibility, which is simpler than considering linear response. 
We shall do this by computing one loop corrections to the classical result in the bulk. 
The flavour of the computation is identical to that for free fields in section [2| The crucial 
difference is that the one loop contribution is to be computed in the curved black hole 
background of the previous subsection, which is 3-1-1 dimensional, as opposed to the 2-1-1 
dimensions of the (strongly coupled) field theory. 

There are several different sources of 1/iV corrections to the free energy. It is helpful 
to identify those most likely to be related to the quantum oscillation structure we are 
seeking. The most universal one loop corrections to the free energy are those coming 
2^ 



We did not specify the connection between the normalisation of the action ( 26 1 and some dual field 
theoretical quantity A''. In general one expects that i^/ft^ scales like A*' to a positive power. That the 
coefflcient of the classical action is large is precisely what allows the bulk side of the AdS/CFT correspondence 
to be treated classically. 
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from the graviton and Maxwell field in (26). These will likely not lead to Landau-level 



related structure, however, as both fields are neutral. The same comment applies to higher 



derivative corrections to the classical action (26). Instead we will focus on the contribution 



of an additional charged field, vanishing in the dyonic black hole background, which could 
be bosonic or fermionic. For bosons the action takes the form 



while for fermions 



(fx^ 



'4)T ■[d+ luJab^ -iqA]il) + mipil) 



(37) 



(38) 



where u^ab is the spin connection. Roman letters denote tangent space indices. 

There are one-loop contributions to the free energy from fluctuations of the scalar and 
fermionic fields: 



9.= VLq + VLb + ^f = ^0 + Ttrlog 



Ttr log 



T -D + m 



+ 



(39) 



where is the classical result (34), V = V — iqA and D = d + jOJabT"'^ — iqA. The boson 



and fermion masses in (39) need not be the same, of course. The dots in (39) indicate that 



we are not computing the one loop contribution from the neutral fields A and g. While the 
classical contribution JIq will scale as some positive power of A^, the one loop logarithms in 



(39) are order one. This is the sense in which we are computing a '1/A^' effect. 



3.3 Determinants in black hole backgrounds and quasinormal modes 

In order to cleanly extract possible T = non-analyticities in the one loop determinants 



(39), we would like to obtain an expression analogous to (13) in the free field case. To do 



this, we must first write down the eigenvalue equation for bosons 



t/^(l) + m^4> = X4>. 



and for fermions 

r • Dijj + mip = Xtp . 
The next step is to separate variables. For bosons this is done by writing 

(j) = e-''^"^+'^yX^{x)(t){r) . 



(40) 
(41) 

(42) 



The quantities appearing in this expression are identical to those in section |2j The important 
difference is that there is one more dimension, the bulk radial direction, and hence a new 
function (/>(r). 
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Separation of variables is a little more complicated for spinors in a magnetic field because 
there are several components that couple differently to the field. However, it is straight- 
forwardly achieved following Feynman and Gell-Mann |30] . One introduces the auxiliary 
spinor x defined by 

^ = (r.£) + A-m)x, (43) 
which is found to satisfy the second order equation 

- b\ + ^-Rx + ^F^^r^^x + (m - \)\ = . (44) 

This second order equation can now be separated exactly as in the bosonic case 

X = e-'^-^+'^yX,(x)x(r) . (45) 



Every x satisfying the second order equation (44) gives an eigenspinor of the original 



Dirac operator (41 ). The solutions will be double counted, because of the extra derivatives. 
However, the matrix commutes with the second order operator. Therefore by imposing, 
say, T^x = X oii6 obtains the correct eigenfunctions without double counting. 

As in section |2] above we will want to analytically continue cOn into the complex plane. 
Setting z = iLOn and substituting the separation of variables ansatze into the eigenvalue 
equations, we obtain 'reduced' equations for each mode 

MB{z,e)(j) = X{z, £)(!}, MpizJ)^ = X{z,i)^. (46) 

These differential equations for (j)(r) and ip^r), when viewed as eigenvalue problems, will 



provide a connection between A and z, similar to (15) in the free case. As previously, the 



k momentum can be eliminated from the equations and only leads to the Landau level 



degeneracy. The important difference between (46) and (15) is that the M^/p are now 
differential operators in the radial direction, so we do not have an algebraic expression for 
A(z,£). 

Mimicking the free theory procedure, the idea now is to express the determinants as sums 
over specific complex frequencies z*(£) that lead to zero modes; A(z^(^),£) = solutions 



of (46). Because M^/p are differential operators we expect to find infinitely many such 
frequencies. For the operators to be well defined, we need to specify the boundary condition 
of (j){r) and ip{r) near the horizon at r = r+ and near the boundary r = 0. The subtler 
boundary condition is at the horizon. The general radial behavior near the horizon is found 
to be 

iz 

0, ~ (r - r+)" + • • • , with a = ±— . (47) 
14 



In computing the Euclidean determinant directly as a sum over eigenvalues, regularity at 
the Euclidean horizon requires taking 

a = • (48) 

This then shows that once we have defined the boundary condition for M^jpi^z^t) on the 
imaginary z = axis, the positive and negative values of LOn will have different analytic 
continuations into the complex z plane. It is important to treat this point carefully in 
deriving the formula we present shortly. 



At general complex z the two boundary conditions in (47) can be called ingoing (the 
minus sign) and outgoing (the positive sign). This corresponds to whether the corresponding 
Lorentzian signature solutions have flux going into the future horizon of the black hole, or 
coming out of the past horizon. On shell modes, with A(zy,(-£),^) = 0, satisfying ingoing 
boundary conditions at the horizon are called quasinormal modes. 

The quasinormal frequencies of a wave equation in a black hole spacetime are poles 
in the corresponding retarded Green's function in the black hole background. To see this 
explicitly it is useful to consider the trace of the inverse of our operators M^jp^ which we 
will denote collectively as M. Starting on the imaginary axis we have 

tr— — -= / G{iuJn,i,r,r)dr, (49) 



M{iuJn,i) Jo 

where the Euclidean Green's function satisfies 



M(iu;„, £) G{iuJn, r, r') = r^6ir, r') . (50) 



The expression (49) follows directly from the usual representation of the Green's function as 



a sum over eigenfunctions. The boundary condition for the Green's function at the horizon 



is (47) together with (48). 



Now consider the analytic continuation of this trace to general complex z = iuJn, where 



we analytically continue (48) from the upper imaginary axis. That is, we take the minus 
sign (ingoing) boundary condition in (47). Denote this object by tr_-p^j-^. In general one 



needs to perform the integral in (49) before analytically continuing. It is clear that the 

poles in this analytically continued Green's function with ingoing boundary conditions at 

the horizon are given by precisely the quasinormal frequencies of the black hole, as this is 

when M{z, i) has a zero eigenvalue. See e.g. [3T] for a more detailed discussionj^ As usual, 

^The quasinormal modes also give the poles of the retarded Green's function of the operator dual to the 
bulk field in the dual field theory [32] ■ The field theory Green's function is essentially given by the behaviour 
of our bulk Green's function near the boundary at r = [35] . 
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continuing the Euclidean Green's function from the upper imaginary axis gives the retarded 
Green's function. If the black hole is stable against linearised perturbations (as they will 
be in the cases we study below) then these poles are necessarily in the lower half plane. 
Furthermore, at finite temperature, the quasinormal modes give isolated poles. 



The conclusion of the previous paragraph is that the nonanalyticities of tr_ 



are 



M{z,e) 

isolated poles in the lower half z plane. We could have instead analytically continued 
the Euclidean Green's function from the negative imaginary axis. Denote this object by 



tr. 



1 

M{z,£) 



. The + boundary condition at the horizon corresponds to outgoing modes. This 



necessarily leads to the advanced Green's function, with poles in the upper half plane. In 
fact 

' ^ (51) 



tr. 



M{z,£) 



tr_ 



1 



M{z, £) 



This relation follows from taking the complex conjugate of (46) and (47). In the following 



we will express our results in terms of the poles of the retarded Green's function (the 
quasinormal modes), as these are more physical for most purposes. If we wish we can always 



obtain the poles of the advanced Green's function from (|51|). 
In the paper 



we derive the following formulae expressing the one loop contributions 
to the action coming from bosonic and fermionic determinants as a sum over the quasinor- 
mal modes of the operators Mb and Mp respectively. The reader may also find appendix [B] 
useful, in which we derive an analogous formula for the simple case of a single damped har- 
monic oscillator. No assumption is made about the quasinormal modes forming a complete 
basis. For bosons 



\qB\AT 
2tt 



e 2*(£) 



4tt^T 



2ttT 



+ Loc . 



(52) 



For fermions we obtain 




(53) 



There difference between bosons and fermions is due to the different thermal frequencies 
([7| and (|8|. In both of these two expressions, the Loc term refers to a 'local' contribution 
to the one loop effective action for the metric and Maxwell fields induced by integrating 
out the charged bosons and fermions. We will discuss these terms a little more below, 
they will not contribute to the various interesting effects we are looking for. Finally, we 



should note that while we have written (52) and (53) in a way adapted to Landau levels 
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and magnetic fields, the representation of determinants in black hole backgrounds as sums 



over quasinormal modes is much more general |28J. The formulae (52) and (53) will be the 
strong coupling analogues of equation ([Ts]) . 

Generally the sums over £ and z^,{£) in (52) and (53) do not converge. These are 
high frequency divergences that should be renormalised, for instance using zeta function 
regularisation. Sometimes to control the asymptotic behavior it is useful to take a step 
back from the above expressions and reintroduce a sum over the thermal frequencies: 

2 



B 



\qB\AT 
2tt 



EE E'"^ 

t z^e) \n>0 



n + 



iz.{£) 



2-kT 



log 



2ttT 



+ Loc . 



(54) 



An entirely analogous expression exists for fermions. This formula is related to the result 

T{z) 



(52) using the following identity from zeta function regularisation: 

1 



^ log(n + z) 



n=0 



d 



ds ^ in + zY 

n=0 



log 



27r 



(55) 



Yet another expression for the determinant is in a spectral representation form. This is 



derived from (54) using contour integration and the fact that the Zi,{£) are all in the lower 
half plane. For bosons we have 

""^ dn 1 



\qB\A^^ 



27r 



Im log (z+ 



n) + Loc . 



(56) 



We will not develop this expression further. A rigorous treatment would need to address 
the validity of closing the contour and the divergences of (56) at 0, = and 0, = —oo. 
Regularity at i7 = may impose constraints on the quasinormal modes. In appendix [B] we 
show how this works for the case of a single damped harmonic oscillator. There is again an 
analogous integral expression for fermions. 

Finally, we should say a few words about the 'local' contribution Loc. Essentially Loc 
contains the local UV counterterms as well as terms that ensure the correct large mass 



behavior in (52) and (53). Equality of the right and left hand sides of these formulae at 



large mass, including order one terms, requires |28j (for bosons, say) 



Loc = < Ttr log 



-V^ + m^j + ^ log 



47r2T 



(57) 



Where |^>o means that we should only keep the terms which remain nonzero in the limit 
A — > oo (the 'nonpolar' terms). Here A determines the scaling of the field near the bound- 
ary, and is related to the mass by the standard AdS/CFT formulae 

3 



A(A - 3) = L^m^ (bosons) , 



A = — h Lm (fermions, m > —h) 



(58) 
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To shorten the expression (57) we have written to include the sum over the Landau 



levels and their degeneracy. The reason that A appears in (57) is that this is the quantity 



that determines the asymptotic boundary conditions. The proof in [28j of the central 
formulae (52) and (53) uses analyticity arguments in A rather than m^. 



The first term in (57) is closely related to the large mass limit of a determinant of the 



form Laplacian plus mass squared. It is well known, see e.g. [M] for a review, that the only 
terms that survive the large mass expansion of such a determinant are given by integrals of 
local curvatures of the background metric and Maxwell fields. Therefore, the effect of this 



first term is to renormalise the Einstein-Maxwell action (26), including the generation of 



higher curvature terms. These terms are blind to Landau levels and therefore will not lead 
to nonanalytic physics as a function of the magnetic field. 



The second term in (57) could likely be computed in principle by using WKB methods 



to obtain the quasinormal frequencies to the first few leading orders in a 1/A expansion, 
perhaps along the lines of [33] . These WKB computations would not be expected to detect 
nonanalyticities of the sort we will describe shortly, which occur at low or zero frequencies. 
In the following we will therefore generally ignore the Loc contribution to the determinant. 



3.4 Zero temperature nonanalyticities 



The zero temperature limit of (52) and (53) is especially simple. As our theory is scale 
invariant, only the ratios -B//^^ and T / fi are meaningful. Let us work at fixed B / ^'^ and 
take the limit T//i 0. How do the quasinormal poles behave in this limit? The two 
possibilities for a given quasinormal frequency are firstly that ^ 0, for instance if 
~ T, and secondly that remains finite, which requires that ~ //. We will see 
explicitly in section |4] below that both possibilities occur. The quasinormal modes that go 
to zero with temperature coalesce and form a branch cut at zero temperature. 



Formally taking the low temperature limit of (52) or (53) gives 



lim ^BIF 



\qB\A 
2-K 



EE 



1. 



Im 



^.(^)lo^ 



2'KT 



+ 



(59) 



In this expression the logarithmic branch cut must be taken along the positive imaginary z 



axis. This is determined by the singularities of the gamma functions in (52 ) and (53) which 



are along the positive imaginary z axis. This zero temperature limit is discussed for the 
damped harmonic oscillator in appendix [B| 

The sum in (59) will only get finite contributions from modes that scale as ~ ^ at 



low temperatures. Frequencies that go to zero with T give a vanishing contribution, as is 
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already discernable in (52) and (53). However, the finite contribution can come from either 
isolated poles or those coalescing to give a branch cut: even though the coalescing poles 
eventually go to zero with T, at any finite T there will be coalescing poles with Zi, ~ For 



the poles forming a branch cut, the sum in (59) will become an integral. 



In general the low temperature sum (59) is still difficult to perform. One difficulty are 
the UV divergences in the sums. We will present in section |4] below some WKB results for 
the large frequency quasinormal modes that are a first step towards a direct evaluation of the 
UV tail of this formula. However, there are specific situations in which the representation 
as a sum of quasinormal modes becomes extremely useful. This is when a particular mode 
or set of modes undergoes nonanalytic motion as a function of a parameter such as Bj 
Derivatives with respect to this parameter will then pick out the contribution of these 
particular modes as dominating over the others. Using results from we will shortly 



perform the sum (59) exactly over a set of poles close to the real frequency axis that 
undergo nonanalytic motion as a function of the magnetic field. 

It was shown in [11 that quasinormal frequencies of charged fermions, i.e. A = 



solutions to the Dirac equation (41) with ingoing boundary conditions at the horizon, can 
undergo nonanalytic motion as a function of spatial momentum k. Specifically, if the charge 
of the fermion is big enough compared to its mass, 3 rr?!? < g^7^, then there exists a critical 
momentum k = kp at which a quasinormal mode bounces off the real frequency axis at 
z = 0. This leads to a low energy peak in the spectral function of the dual field theory 
fermionic operator near to a particular finite momentum. At T = the peaks becomes 
a delta function. The momentum kp was therefore identified as the 'Fermi momentum' 
indicative of an underling strongly coupled Fermi surface. 

The results from [11], at finite momentum but zero magnetic field, can be adapted to 
our context as follows]^ We can note that the magnetic field B appears in the 'second 
order Dirac equation' ( |44[ ) in two ways. Firstly it appears as just B in the metric function 
/(r) and in the spin-magnetic 'Zeeman' interation FabT"'^- Secondly, it appears as £B, i.e. 
multiplied by the Landau level, in the gauge covariant kinetic term. If we take the limit 
S — > with 2£\qB\ = k'^ fixed then we loose the first terms while retaining the kinetic 



term. As in the free field case discussed around ( 18 ) above, this limit reproduces precisely 
the B = and finite momentum k equation studied in [11 . We can therefore directly 
use results from that paper, with the pole now bouncing off the real axis at 2i\qB\ = kp. 
The i? — > with IB fixed limit is not essential to use results from [Tl]. Keeping B finite 



* Fermionic quasinormal modes in a magnetic field were recently studied in [361 [3] 
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introduces some smooth B dependence into the various 'constants' that appear in this and 
the following sections. 



3.5 Summing low temperature poles: Quantum oscillations 

Before taking the strict zero temperature limit, it is useful to look at the pole motion 
at finite but low temperature. At frequencies and temperatures that are small compared 
to the chemical potential, z,T <^ fi, it is possible to solve the Dirac equation explicitly, 
see appendix D4 of [11]. Using the observation of the previous section we may translate 
the expressions from that paper into results for the quasinormal frequencies with a finite 
magnetic field in the limit B ^ with iB fixed. As we noted, this limit is not essential but 
cleanly extracts the nonanalytic behavior. 

A crucial parameter in the discussion of [TP is u. This quantity controls the low energy 
{uj <^ fi) scaling dimension of the dual fermionic operator in the strongly coupled field 
theory. This scaling dimension is related to the charge and mass of the field by 




The Fermi momentum in units of the chemical potential, 'ykp/n, also depends on m and 
q. This dependence must be determined by numerically solving the Dirac equation in the 
Reissner-Nordstrom black hole background. A plot of z/ as a function of m and q may be 
found in figure 6 of [11] . It can be shown that u is always real. 

We will assume for concreteness that u < ^ {a similar discussion will hold for the case 
u > In this case, the quasinormal frequencies z^, in the low temperature and small 
frequency regime were found to be given by 

Hz.) = , (61) 

where 

. (62) 



We have rearranged the expression appearing in [TT] because it will be important that J-{z) 



has zeros but no poles. In (62) we have introduced 

= ^/2^\qB\ - kp , (63) 
which is a measure of the magnetic field and can be either positive or negative. The 



constants h and 6 in ( 62 1 are determined in jllj in terms of the charge and mass of the 
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fermionic field (numerically in the case of h). It will be sufficient for our purposes to take 
them to be order one in units of the chemical potential. The value of 6 is constrained to lie 
in the range 

< 6* < 7r(l - 2z^) , (64) 
which guarantees that the poles are in the lower half plane for both signs of k±. 



The equation (61) will clearly lead to quasinormal frequencies of the form 



(65) 



for some sequence of functions J^^^\ It is straightforward to solve (61) numerically and 



obtain the motion of the quasinormal poles as a function of k± In figure [s] we show the 
low temperature motion of the poles closest to the real axis as k± is varied through zero, 
for a particular choice of numerical values of the parameters involved. 

0.0 r 




0.0 

Re(z//i) 



Figure 3: Motion of the quasinormal frequencies closest to the real axis as k±/ fi is varied 



from —1 to +1, according to (61). The temperature is T = 0.005/i. The other constants 



are taken to have values g = 1, 7 = y/l2, u = 1/3,9 = 7r/6,h = fx^^^. 

In figure [3] we see several interesting effects. Firstly we can see the advertised pole that 
moves up to and then sharply bounces off the real axis. The bounce has been smoothed 
out at finite temperature. Secondly, there are poles coalescing to form a zero temperature 
branch cut. These poles show a nontrivial circular motion as a function of k±. We now 



^The authors of [TT] have considered this problem in detail. We thank John McGreevy for drawing our 
attention to their appendix D4 and for sharing unpublished results on the motion of quasinormal poles in 
this low frequency regime. 
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need to compute the magnitude of these effects on the magnetic susceptibiHty as k± goes 
through zero. 

Figure |4] shows the contribution of these lowest few quasinormal poles to the magnetic 



susceptibility as a function of k± . These are computed using our formula ( 53 ) and strictly 



speaking we plot the quantity x, see (68) below, which is closely related to the susceptibility. 



The darker line in the first plot is the contribution of the 'T = 0' pole that bounces off 
the real axis in [3| The figure also shows the total susceptibility arising from the sum of 
the contributions of the lowest fifty modes. As anticipated, there is a strong feature in the 
response around = 0. We can also discern other features in the individual responses 
of the modes. Somewhat magically, the motion of the 'branch cut' poles is choreographed 
to precisely cancel out these extra features. The contribution of these other poles are the 
lighter lines in the first plot, while the second (right hand) plot shows the total response due 
to the lowest fifty poles. In the second plot only a single feature remains in the magnetic 
response. The cancellation between oscillations may perhaps be thought of as analogous to 
a Fourier transformation, in which sums of oscillations can cancel to give simple functions. 





Figure 4: Left: Contributions of the lowest few quasinormal modes to the magnetic sus- 



ceptibility, according to (53), as a function of kj_/fi. The darker line is the pole nearest 
the real axis. Right: the total magnetic susceptibility due to the lowest fifty modes. The 
temperature is T = 0.005^. The constants have the same values as in figure |3j The vertical 



axis is proportional to x of ( 68 ) . 



The peak seen in figure gwih occur whenever 2£\qB\ = kp. Thus the peaks are periodic 
m 1/B with period 2TTq/AF, as expected for quantum oscillations due to a Fermi surface. 
We will shortly make this statement sharper by going to the zero temperature limit. 

The right hand plot in figure l4l the sum of the lowest fifty poles, only makes sense if 
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the series being summed is convergent. We are interested in the magnetic susceptibihty, 



X = —5^0, with given by (59) and the quasinormal modes Zj, given by (61 ). To determine 



convergence of this second derivative of the sum we need to know the dependence of on 
the magnetic field B at large values of z*. Let us focus on B close to the critical value Bi 



at which the ^th oscillation occurs: B = B^ + SB = kp/2£q + dB. Then from (63) 



k^ = ^SB + ... . 



(66) 



From (61 ), by expanding the right hand gamma function in (62 ) in the vicinity of a negative 



integer, we find that at large Zi, and for these small values of k_\_: 

dz^ 



dB 



-2u 



(67) 



We can now see by differentiating (59) and using (67) that while the sum over Zi,{i) in 
X = —dl^U, is UV divergent, this is only due to derivatives acting on the 'trivial' overall 



factor of B in (59). The divergent factor can be removed by considering for instance 



—d^Vl + 2dB{S^B~^)^ which does lead to a convergent sum. As we expect the term with 
most derivatives of O with respect to B to capture the strongest nonanalyticities, extra terms 
depending on single derivatives of O should not be important for k±_ ~ 0. Alternatively we 



can define, using (66), 



X 



ki dki ' 



(68) 



where is the ith component of 17 = B^^Qg. This leads to a convergent sum over -z*(^) 
and is again equivalent to x up to first derivatives of 0,. In particular, we expect x ^ X 
very low temperatures and A;_l ~ 0. The finite quantity x been used as the vertical axis 
of figure |4| 

Having obtained convergent sums over Zi,{l), the sum over i itself is still not convergent. 
We shall not be concerned with this divergence, as we are considering low temperature 
nonanalyticities that occur for each i individually at different values of the magnetic field. 
These nonanalyticities are not sensitive to the large i UV divergences, analogously to the 
free field case of section|2j In order to exhibit the quantum oscillations at high temperatures, 
one will likely have to perform the sum over i. 



3.6 The zero temperature limit 

In this section we will obtain the susceptibility exactly at zero temperature and for kj 



0. 



To this end, the sum over quasinormal modes in (59) is helpfully rewritten as an integral 
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along the real frequency axis 

\qB\A 



2tt 



J^-Ini ^ 



vr 



2Tri 



- iz T'{z) 

_^ ^ 27rr ^(z) 



(69) 



This expression follows from contour integration and the fact that given in (62), has 



zeros at the quasinormal modes and no poles. The expression (69) is somewhat formal, 



but we now take two derivatives with respect to B to obtain a convergent expression as 
described at the end of the previous section. 



The zero temperature limit of the susceptibility is obtained by differentiating ( 69 ) . Using 



( 68 ) for the susceptibility (recall that x ~ X ™ this regime) : 



\3B\AS. 

27r ^ 27r2/c2 

f F 



Re 



{k± - he^^z'^'^ f 



log — "^dz . 



(70) 



In the integrand it is important to take the branch cut due to the powers z"^^ to run down the 
negative imaginary axis, this is required by the coalescence of poles of the gamma functions 



in ( 62 ) . As we noted previously, the logarithmic branch cut must run along the positive 



imaginary axis. The integral can be performed exactly to yield 
^ \qB\A{2v-l) q^ =ir. A _ o2 



sm 



X 



27r 



4z^2 klh^/^u sin ^ ^ (_A;^)2-i/2^ 



(A;± < 0) 



(71) 



For A;_|_ > one replaces k\ 



-k I and sin 



2v 



sm 



2v • 



The integral is only convergent 



if 4 < z/. This extra condition is required to be able to close the contour in the lower half 



plane in the derivation of (70) and is also the condition for the power of k± appearing in 



(71 ) to be negative. For smaller values of v one needs to differentiate the free energy more 



times to obtain a convergent integral and furthermore a divergent dependence on k±. There 



is no temperature dependence in (71). Technically this occurs because the integral in (70) 



vanishes if the logarithmic term is not included. This shows that the log T in ( |70[ ) does 
not lead to a logarithmic divergence in the susceptibility at low temperatures]^ The zero 
temperature result ( 71 ) is plotted in figure [5} 



Schematically, ( 71 ) can be written as 



X 



lim ^ „o 
T^o dB^ 



+ \qB\A 



2i\qB\ 



kl 



-2+1/21/ 



(72) 



^Vanishing of the integral without the logarithm also indicates that J]]^ z*, suitably regularised, is an 
analytic expression even though individual poles undergo nonanalytic motion. We suspect this may be a 
general phenomenon. Thus Ylz z*logz* is needed to extract the nonanalytic dependence on B. To obtain 
the correct answer one must sum all the poles near the real axis, it is not sufficient to focus on a single pole. 
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Figure 5: The magnetic susceptibility at T = as a function of k±/fi. The constants have 
the same values as in figure [Sj The vertical axis is proportional to x of ( 68 ) . 



The sign is important and physical. The divergences in the susceptibility at 2i\qB\ = kp 
are seen to be positive, with opposite sign to the delta functions appearing for free fermions 



in (21). The sign follows from the observation that sin J^/ sin < in the region we 



are studying: 1 < u < I and (64). The sum over £ can be performed in (72) in terms of 



generalised zeta functions. The formula ( 72 ) is analogous to the result (21 ) for free fermions. 



Once again, it indicates the existence of oscillations in the magnetic susceptibility with 
period 



A 



(73) 



1 \ _ 27rg 

b) ~a:^- 

As well as the sign with which the susceptibility diverges, another importance difference 



with respect to the free fermion result (21 ) is that the nonanalyticity is softer in the strongly 



coupled theory. Rather than the delta functions of the free theory (21 ) we find the absolute 



value of a (generally non-integer) power in ( 72 ) . The power is determined by the low energy 



scaling dimension u in (60). Our computation is valid for 4 < u < I. These inequalities 



are satisfied for a range of values of q and m, including for instance m = with jq = 1, see 
figure 6 of [TT] . 

To restate the main results of the last few sections 

• For a range of values of the mass m and charge q of the bulk fermion there is a 
quasinormal pole which (at T = 0) nonanalytically bounces off the origin of the real 
frequency axis at 2i\qB\ = kp. 
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For a certain range of m and q, corresponding to | < < |, these bounces produce 
periodic in 1/B divergences in the one loop magnetic susceptibihty. The periodicity 



is given by ( 73 ) and the strength of the divergence by ( 72 ) . 



This behavior would seem to be aptly characterised as a strong coupling manifestation of de 
Haas-van Alphen oscillations at low temperatures. Thus our results simultaneously support 
the characterisation of kp as a Fermi momentum in [HI |9l [101 [H] ^-Iso indicate qualitative 
differences between de Haas-van Alphen oscillations at weak and strong coupling. 

We now turn to a numerical study of the quasinormal modes of charged bosons and 
show that the modes can have an interesting magnetic field dependence in that case also. 
Quantum oscillations from stable charged bosons would be a novel effect. 

4 Charged scalar quasinormal modes of dyonic black holes 
4.1 Equations for bosons 

In the previous section we found analytic results for quantum oscillations due to fermionic 
quasinormal modes. This was possible because the relevant nonanalytic motion of the 
quasinormal mode occurred close to zero frequency, ~ 0, as the mode bounced off the 
real frequency axis. The Dirac equation in the zero temperature AdS-Reissner-Nordstrom 
black hole was solved analytically at small frequencies in [11 J. 

For bosons, in contrast, if a quasinormal mode moves towards z^, ~ it typically indicates 
the onset of a superconducting instability jH E]- Rather than bounce back into the 
lower half frequency plane, the mode continues up into the upper half plane causing an 
instability and the condensation of the bosonic field. Therefore, if we wish to look for 
possible nonanalytic motion of bosonic quasinormal modes, without going through a phase 
transition, we will need to look away from small frequencies. To this end we will study the 
bosonic quasinormal modes numerically. The hunt for analogues of quantum oscillations 
leads us to look for special values of Ki ~ £B. Near the end of this section we will also look 
at the magnetic susceptibility at the onset of superconductivity. 



Recall that after separating variables as in (42), and analytically continuing z = iuJn, 
the eigenvalue equation became 

Mb{z,£)^ = \(t). (74) 
The 'reduced' operator is found to take the form 

26 



Recall that / was given in equation (30). In Appendix A we put this eigenvalue equation 
in Schrodinger form. We are looking for the quasinormal modes of the operator Mb- These 
are A = eigenmodes of Mb satisfying ingoing boundary conditions at the horizon and 
normalisable boundary conditions at infinity. We have already noted that ingoing boundary 



conditions at the horizon (r = r-|_) corresponds to taking the minus sign in (47). Near the 



asymptotic boundary of the spacetime (r = 0) the general behavior of A = modes is 



^ = + • • • , with a = ^ ± y ^ + L'^m'^ . (76) 

Normalisability at infinity generally requires taking the faster of the two falloffs. In this 
paper we will ignore the possible ambiguities that arise for masses sufficiently close to the 
Breitenlohner-Freedman bound {jn?BF ~ ~9/4L^) and simply impose the faster falloff at 
the boundary. For a discussion of determinants in AdS with m^p < rn?' < rrigp + 1 see e.g. 



4.2 The matrix method for quasinormal modes 

The quasinormal modes of a black hole are complex frequencies z^{£) such that there are 
solutions (j) satisfying 

Mb(z.(^),^)0 = O, (77) 
together with ingoing boundary conditions at the horizon and normalisability at infinity. 



These are the frequencies that contribute to our sum ( 52 ) . 

The technical challenge we face is to find the quasinormal modes of a charged scalar field 
in planar dyonic Reissner-Nordstrom-AdS black holes. In particular, we will be interested in 
low and zero temperatures. While there is an immense literature on quasinormal modes, to 
our knowledge this particular problem has not been addressed. The most relevant references 
are collected in section 6.2 of the review [40J. It was noted in [Jl] that sometimes quasi- 
normal modes in AdS are easier to find than in asymptotically flat spacetimes, because the 
AdS conformal boundary gives a regular singular point in the relevant differential equation 
rather than an irregular singular point. Unfortunately, the techniques of |3T] will not work 



for us because at low temperatures there are singular points in the differential equation ( 77 ) 
that become arbitrary close to the horizon (r = r-(_) and make a Taylor series expansion at 
the horizon useless. At strictly zero temperature the horizon becomes an essential singular 
point and a series expansion there has zero radius of convergence. 

A useful discussion of asymptotically flat zero temperature Reissner-Nordstrom quasi- 
normal modes can be found in [42j. The authors of that paper noted that, after a change 
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of variables to bring infinity to a finite radial coordinate, then a Taylor series expansion 
at the midpoint between the horizon and infinity had a radius of convergence that reached 



both the horizon and infinity. The same property holds for our equation (77): the Taylor 
series about the midpoint rmid = has radius of convergence and therefore reaches 
both the horizon at r = r-|_ and the boundary r = 0. This is true for all values of various 
parameters in the equation, including the zero temperature limit. 

While |32] were then able to find the quasinormal modes by reducing a 5-term recurrence 
relation for the Taylor series about the midpoint to two 3-term recurrence relations and then 
using a continued fraction method due to Leaver |33j, our case is more complicated. A Taylor 



series expansion of (77) about r = leads to a 9-term recurrence relation. Fortunately 
|43j also presented a method for dealing with arbitrary length recurrence relations. We will 
now review the algorithm. 

1. Expand (/> in a series expansion about the midpoint, having first factored out the 
desired leading (singular) behavior at the horizon and infinity. Thus at finite temper- 
ature 

1 ^ 



n=0 

while at zero temperature 

N 



Jzrl/6(r+-r) J-i(4:Z-3q^,)r+/36^^(3+^/9+4LW) ^ ^^^^ _ l^^^n ^ ^jg^^ 



n=0 
3 , o /„ \4 



Note that at zero temperature /(r) = 1 — 4 {r/r^) + 3 {r/rj^) 



2. Plug the relevant series expansion into the differential equation (77) and expand. 
Collecting in powers of r — ^r-|_ gives A^+l linear relations between the A^+l coefficients 
{an}- Write these as a matrix equation: 

N 

Y,Amn{z)an = ^. (80) 
n=0 

Because the recurrence relation between the {an} involves nine terms in general, A 
will have nonzero entries only along a diagonal band of width nine. 

3. The quasinormal modes are given by the zeros of the determinant of the matrix A 

det A{z^) = . (81) 

Given that the matrix A is fairly sparse, this determinant can be numerically computed 
quickly and robustly using, for instance, Mathematica. 
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As an illustration and to introduce concepts we first present the results of this method 
for a neutral scalar field {q = 0) with no magnetic field background {B = 0) at low and 
zero temperature. For the moment we will make the choice of mass m? = 0; with this mass, 
neutral scalar fields are stable all the way down to zero temperature [?]• Furthermore, for 
concreteness we will take 7 = 1 throughout the remainder of this paper. For some values 
of 7 obtained via Freund-Rubin compactifications of M theory, see [7] . 

Figure |6] shows the quasinormal modes closest to the real axis for a small (left) and 
zero (right) temperature. The structure we are about to describe was to a large extent 
previously noted in |44| H5] . One sees clearly that at finite low temperature there are 
two distinct types of quasinormal modes. Along the negative imaginary axis we have a 
sequence of closely spaced modes, while on each side of the negative imaginary axis there is 
another sequence of modes descending diagonally. Because we are at low temperatures, it is 
natural to think of the closely spaced modes as having positions dominantly determined by 
T whereas the off-axis modes are more sensitive to /i. This statement can be made precise 
by varying the temperature, but we shall not go into detail here. The symmetry of the plot 
under z 1— > — z follows from the corresponding transformation of the differential equation 



(77) and ingoing boundary conditions (47) when q = 0. 

In the zero temperature limit, T//x — > 0, we should expect the modes along the negative 
real axis to bunch together and possibly form a branch cut. The right hand plot in figure 
[6j showing the zero temperature quasinormal modes, supports this picture. The fact that 
discrete poles are still visible in the plot is an artifact of truncating the differential equation 
to a finite matrix equation (with rank A^-|- 1) in step 2 of the algorithm we presented above. 
As N is increased one can check that the modes in the figure move up the imaginary axis, 
becoming arbitrarily clumped as ^ 00. The modes away from the imaginary axis remain 
fixed, consistent with the notion that their spacing is dominantly set by fi. 

For the case of a neutral scalar it is straightforward to argue analytically that there is 
a branch cut in the retarded Green's function at zero temperature (while the quasinormal 
modes correspond to poles in the Green's function) . Making the change of variables 

djT 

</) = rM/ , J = ds, (82) 



the differential equation ( 77 ) with m =g = OatT = becomes the Schrodinger equation 

d^^^mnjn_f^r))^ = u^^, (83) 



where we think of r as r{s). By integrating (82) one easily finds that s +00 as r — > r+. 



Furthermore, the potential in the Schrodinger equation (83) has the leading near-horizon 
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Figure 6: Left: quasinormal modes at temperature T//i = 0.075. Right: quasinormal 
modes at zero temperature T = 0. The plots show the lower half z/fi frequency plane, 
and bright spots denote quasinormal modes. Both plots have q = 0, m? =0, 7 = 1 and 
no magnetic field. The plots were generated as density plots of |det ^'(2:)/ det A(2;)|. The 
finite temperature plot has N = 200 while the zero temperature plot has = 500. The 
discreteness of the poles on the imaginary axis at T = is an artifact of finite A^. As 
N ^ 00 the poles coalesce and form a branch cut. 



The fact that the Schrodinger potential has a power law rather than exponential falloff near 
the horizon is characteristic of extremal rather than finite temperature horizons. When 
a Schrodinger equation has an asymptotic region in which the potential has power law 
falloff, one generically expects that the late time evolution of the wavefunction at some 
fixed position will be dominated by scattering events in which an excitation travels a long 
distance into the asymptotic region and is reflected back. This occurs at a rate proportional 
to the asymptotic potential and hence leads to a power law decay in time (see e.g. |l6l|31] for 
more precise arguments). Fourier transforming to frequency space, the power law tail leads 
to a branch cut running along the negative imaginary axis from z = 0. This is the branch 
cut we are seeing in figure |6] at T = 0. In contrast, an exponential falloff near the horizon can 
be shown to imply that there cannot exist branch cuts emanating from z = |47| US] . The 
presence of a late time power law tail in extreme Reissner-Nordstrom-AdS was previously 
shown in ^45^, in agreement with the observations we have just made. Furthermore 
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(84) 
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have recently exhibited this branch cut by exphcitly solving the Schrodinger equation near 
z = 0. 



4.3 Charged scalars and magnetic field dependence 
4.3.1 The scales involved: T,B,n 

We now wish to determine the behavior of the quasinormal modes as a function of the 
magnetic field in the low temperature limit. There are three scales characterising the 
black hole and dual field theory: T, /i and B. The remaining dimensionful quantity is the 
quasinormal mode frequency z^. Because there are no other scales, the underlying strongly 
coupled theory is a CFT, only the ratios of these dimensionful quantities can be physical. 
We will implement this freedom as follows. 

Firstly, note that the following rescalings completely eliminate r_|_ and 7 from our dif- 



ferential equation ( 77 1 



r = — , z = zr+, T = Tr+ , /x = , B = , q = q-f . (85) 

r+ 7 7 



The quantities {T, are dimensionless and satisfy the constraint, from (33) 



4ttT = 3- fi^ - B\ (86) 

We will use this constraint to eliminate, for instance, fi. We will then find the dimensionless 
quasinormal modes 3,S Oil function of T and B using the algorithm above. Finally, the 
physical quasinormal frequency in units of the chemical potential is obtained by 

- = 4 = , ^ . (87) 

7/^ 7x/3 - 52 - 47rf 

This will be a quasinormal mode at physical temperature and magnetic field 

^= , ^ 4 = ^ ^. (88) 

The upshot of the considerations of the previous paragraph will be quasinormal modes 
of the form 

,ip)^i)=:r(^)^l^^j^^^ (89) 

for some sequence of functions T^^h Our objective is to explore these functions, particularly 



the B dependence. This will then allow us to evaluate some of the terms in the sum (52). 
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4.3.2 Dependence on the charge q of the scalar field 

Before switching on a magnetic field we make a few observations about the q dependence. 
Without loss of generality we will focus on positive charges. Taking negative charge q would 
simply result in a reflection about the imaginary axis: 

q ^ -q <^ ^ -Zi, . (90) 



This follows from taking the complex conjugate of (75 1 together with noting that z —z 
will preserve ingoing boundary conditions at the horizon. 

For a given mass squared there exists a critical charge such that if q is larger than 
the critical value the extremal Reissner-Nordstrom black hole becomes unstable. A precise 
expression for the critical charge as a function of the mass was obtained in |7] , following the 
initial discussions of the instability in [H [5l |6]. New instabilities were discovered in 
we briefly discuss these near the end of this section. The instability indicates the onset of 
a superconducting phase. A finite temperature improves the stability of the black holes, 
although a sufficiently large charge q will always result in an instability. In figure [7] we show 
how the location of the quasinormal mode closest to the real axis moves as a function of 
charge g at a fixed low temperature T/^ = 0.05 and m? = 0. At g = the mode is on 
the imaginary axis, as in figure |6] As the charge is increased the mode moves up towards 
the real axis following an almost semicircular trajectory. At a critical charge qc ~ 4.3 the 
pole crosses into the upper half plane, indicating the onset of a superconducting instability. 
This value for q^ agrees nicely with figure 1 in [7 . 

While the onset of superconductivity is the most dramatic effect that occurs as a function 
of the charge of the scalar field, it is also interesting to see how the higher quasinormal modes 
rearrange themselves. The low lying quasinormal modes at low temperature of a scalar field 
with charges q = 1,2 and 4 are shown in figure [8j The g = quasinormal modes for the 
same mass and temperature were already shown in figure |6j 

In the plots of figure |8] we can see how the mode closest to the real axis moves up towards 
the axis as the charge is increased, as we saw previously in figure [7j The motion of the 
other quasinormal modes is very curious. The line of quasinormal modes that was along 
the negative imaginary axis at zero charge bends increasingly towards the left. Meanwhile, 
the modes that were to the left of these in figure |6] move upwards and to the right, and 
one by one merge with the sequence of modes that were along the imaginary axis. By the 
time the charge is g = 4, the rightmost plot in figure [Sj it is impossible to distinguish any 
more between these two types of low lying modes. On the other hand, the modes that were 
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Figure 7: The location of the quasinormal mode closest to the real axis as a function of 
charge. The temperature is T//i = 0.05 and B = 0. The charge ranges from g = to 
q = 4.3. The mass is = while 7 = 1. At the upper limit of q the mode crosses into the 
upper half plane, indicating the onset of a superconducting instability. 



to the right of the imaginary axis are pushed down and further to the right, eventually 
disappearing from our plot. 

Figure |8] suggests the following interpretation. There are quasinormal modes that for 
charges less than some critical charge q < gcrit. remain of order /i as T ^ 0, while at 
sufficiently large charge q > gcrit. they coalesce with other poles that are forming a branch 
cut as T — > 0. At strictly zero temperature this would presumably correspond to a critical 
charge at which the pole crosses the branch cut into an 'unphysical' sheet. Any given pole 
forming the branch tends to the origin as T — > 0. Therefore one would have the nonanalytic 
in q behavior that 

0{fi) for q < q„it ^^^^ 
for g > Qcrit 

To establish this behavior unambiguously would require higher precision numerics at low 
temperature than we are currently able to perform. 

Varying the charge q corresponds to varying the charge of an operator in the field theory 
and so is not an operation that can be performed within a given theory. However, we will 
now see that the merging effect shown in figure [s] and equation (91) can be undone by a 
magnetic field. This leads to a possible source of (periodic in 1/B) nonanalytic behavior of 
the free energy due to bosons. 
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Figure 8: Quasinormal modes at temperature T/fx = 0.075. Charges from left to right: 
q = 1, q = 2 and q = 4. The plots show the lower half frequency plane z/fi, and bright 
spots denote quasinormal modes. Both plots have m? = 0, 7 = 1 and no magnetic field. 



4.3.3 Dependence on B: possible periodic nonanalyticities from bosons 

We now turn to the dependence of the quasinormal modes on the magnetic field. We will 
focus on a case in which, in the absence of a magnetic field, the merging of poles with a 
branch cut uncovered in figure [8] has occurred. This requires a sufficiently large charge q. 
As we would like to work at zero temperature, we will increase the mass of the field so that 
it remains stable at zero temperature. The choice q = 4:, m? = 10 (hence A = 5) does the 
job [7|. We will also work in the limit of large £ and small B, with BH order one, so that 
we can isolate effects that are potentially periodic in 1/B. 

Figure |9] shows the behavior of the T = quasinormal modes of this charged scalar 
{q = A,m'^ = 10) as a magnetic field is turned on. The top left plot, with a small value of 
iB^ is analogous to the rightmost plot in figure |8j The branch cut bends to the left whereas 
various poles have either merged with the branch cut or moved off to the right of the plotted 
region (as in figure [s]). As the magnetic field is increased, the poles that had merged with 
the branch cut are seen to reappear, moving out to the left. 

Within the current accuracy of our numerics, it is not completely clear whether the poles 
emerge from the cut at specific finite values of iB or whether the pole is distinct from the 
cut at arbitrarily small iB. If the former case is true, than this may lead to non-analyticities 
in the free energy that are periodic in 1/-B, with qcnt. in (91) replaced by a critical A;crit. 
(with k"^ = 2l\qB\). Given that such nonanalyticities are usually associated strictly with 
fermions, it would be very interesting indeed if bosonic operators can lead to these effects 
at strong coupling. We hope that future work will settle this question. 

A possible simple interpretation for why the poles move to the left in figure [9] suggests 
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Figure 9: Emergence of quasinormal modes from a branch cut as a function of magnetic 
field and at zero temperature. All plots have g = 4, 7 = 1, m} = 10 and I = 10. The plots 
show the lower half frequency plane and bright spots denote quasinormal modes. From 
top left kill = y'2£\qB\/fi w 0.283,0.632,0.894,1.26,1.79,2.53. The closely spaced modes 



are a branch cut that has been discretised into poles by the finite value of = 500 in (80). 



itself. Recall that in the limit we are working here, the poles in a magnetic field B at some 
Landau level i are the same as the poles in the absence of a magnetic field and at a spatial 
momentum k"^ = 2i\qB\. At zero chemical potential (/i = 0), relativistic invariance will 
require a branch cut in the T = retarded Green's function emanating from k = z. It may 
be that the vertically spaced poles in the lower rightmost plots of figure |9] are a finite n 
remnant of this branch cut. 



4.4 Magnetic susceptibility near a superconducting instability 

We have already mentioned that if a bosonic quasinormal mode hits the real axis it typically 
indicates the onset of a superconducting instability. We saw an example of this in figure 
([7| as a function of the charge of the scalar field. One expects to see a strong response in 
the susceptibility as the temperature is lowered close to the critical temperature Tc, yet this 



in absent to leading order at large N (i.e. in the classical result of equation (34) which is 
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unaware of the existence of Tc). In this section we will outline a one loop computation of a 
divergent magnetic susceptibility. The computation is essentially identical to the standard 
computation in flat space, see e.g. |39], except that we will use zeta function regularisation. 
The simple steps below form the conceptual outline of a more complete calculation. 

Firstly, we assume that only quasinormal modes close to the origin in frequency space 
are important. This is reasonable because the divergence we find is directly due to the log 
divergence of our formula (52) as z^, — > 0. At least one such mode exists, by assumption. 



because T T^. To obtain the correct susceptibility we must sum over the higher Landau 
levels associated to this mode [49^ . These are also close to the origin in the limit of small 
magnetic field, a limit we will assume in this calculation. In fact, the computation we are 
about to present will require 

5 « (T - n)T, . (92) 

However, we will not sum over the 'excited' quasinormal modes, as these will be at a distance 
at least of order Tc from the origin. Using the fact that T(x) ~l/x + -- -,asx— >0, then 
for the modes of interest, ~ 0, (|52|) reduces to 
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(93) 



Secondly, we assume that the quasinormal modes appearing in (93) can be written as 

z^{£) = ai\qB\{i + i) + Q2(r - T,) + • • • , (94) 

with ai complex constants with negative imaginary parts. The temperature dependence is 
fixed by the requirement that go to zero at T = Tc in the absence of a magnetic field. 
The magnetic field and Landau level dependence at a linearised level at small B is fixed by 



the fact that only the combination B{£+ i) appears in the differential operator (75). The 



offset of 2 is crucial in this computation. Thus the free energy (93) becomes, using zeta 
function regularisation. 
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(95) 
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where as previously (^h{s,x) is the Hurwitz zeta function. 

To leading order in B <^ (T — Tc)Tc this gives a divergent susceptibility 
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(97) 
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If the real part of axja^ is positive, then the divergence is diamagnetic, as we should expect 
for the onset of superconductivity. Numerical computations suggest that this is the case. 
As previously, the divergence indicates the breakdown of perturbation theory and the need 
to resum higher order corrections. This divergence is a second vivid example of physics 
beyond the classical gravity limit, albeit not the focus of the present paper. It is of interest 
to flesh out the computation outlined above. 

A different order of limits is to take T = with a large magnetic field stabilising the 
superconducting instability. One then lowers the magnetic field, finding an instability at 
= Bc2 jSOl |6] . Distinct to the case we have just treated, this would be a zero temperature 
quantum phase transition. Furthermore the zero temperature ground state for B < Bc2 is 



not known in this regime. It is plausible that we can use the same logic as in sections 3.4 



3.6, The relevant (likely £ = 0) quasinormal mode should behave for B just above Bc2 like 

lmz,{0) ^ -\B- B,2f , (98) 



for some positive 6 > 0. The bosonic analogue of (72) will then give 

X ±\qB\A\B - B,2\^-\ (99) 



Interestingly, the sign depends on the value of the angle 6 of section 3.6 which now lies in 
the range 7r(l — 2iy) < 9 < it. For cases in which the divergence is paramagnetic rather than 
diamagnetic, it may be that the state at i? < Bc2 is not superconducting at T = 0. To 
fully address this question one should generalise the results of [11^ for scalars in a magnetic 
field. We cannot simply use the existing results as these translate into statements about 
large £ modes, while the T = finite B instability is likely sensitive to the lowest Landau 
level only, with the higher levels being 'gapped' by the magnetic field. 

4.5 WKB results for charged scalar quasinormal modes 

As a check on our numerics, and also with a view to eventually controlling the UV behavior 



of the sums over quasnormal modes in (52 ) and (53), we have computed the large frequency 
limit of the charged bosonic quasinormal modes, at fixed £. The derivation is given in 
appendix |C] closely following the WKB analysis of neutral massless scalars in the AdS- 
Reissner-Nordstrom background in (51]. The method used for obtaining the asymptotic 
behavior of quasinormal modes was pioneered by [52]. 



The result is that the asymptotic quasinormal modes for the equation ( 75 ) are 



^0 / (r+-r(x))\ 1, / /7r\\ 7r(2 + \/4m2 + 9) , , 
dxU-g/x ^ + " j =n^ + -log(2cos(-jj + ^ ^ (100) 
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with n G N, while xq and x{r) are given by 



--g74y>-(-^). -«^-E74y.o.(i--). (101) 



where Vp are the zeroes of /(r). Thus r(x) in (100) must be found by inverting the second 
of these expressions. The integral in ( 100| ) is along the Stokes line specified in |5lJ. In the 
complex x plane this is simply a straight line connecting and xq. In computing r(x) one 
must be careful to remain on the correct sheet. The 'mirror' poles, z^, —z^ and q — > —q, 
are obtained with a different choice of Stokes lines. 



The result (100) does not depend on the Landau level I and so only depends on B 
through f{r). To capture the Landau level dependence will require a more refined analysis, 
perhaps along the lines of |35j. The complication is that the large £ and large n limits do 
not commute. Ideally we would like a single formula encapsulating both of these limits. 



We have checked that (100) agrees with our numerical results for the sequence of poles 
extending down in the complex plane with a spacing set by /i (the slope agrees very well, 
there is a slight mismatch in the offset which we believe is due to not being able to access the 
large frequency limit with our numerics). A curious feature, however, is that the sequence 
of poles coalescing to form a branch cut, i.e. those that appear more closely spaced at low 
temperatures in figures ([6]) and ([s]), do not appear in this formula. This may suggest that 
the zero temperature branch cut emanating from the origin terminates at some large but 
finite frequency rather than extending asymptotically. We hope this can be determined 
precisely in future work. 



5 Discussion 

In this paper we have argued that one loop effects in the bulk contain information that needs 
to be accessed if the 'applied holography' research program is to successfully disentangle 
the different physical contributions to observables such as the free energy and electrical 
conductivity. We illustrated this point of view by showing that the one loop contribution 
to the magnetic susceptibility exhibits the de Haas-van Alphen quantum oscillations that 
are not manifest at a bulk classical level. In achieving this, we discovered that the oscilla- 
tions at strong coupling differ from their weak coupling counterparts in that the periodic 
low temperature nonanalyticities of the magnetic susceptibility are not delta functions but 
rather power law divergences. 

The key step in our arguments has been the (new to our knowledge) observation that 
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one loop determinants can be expressed as a sum over quasinormal modes. The quasinormal 
frequencies give the poles of retarded Green's functions in the strongly coupled dual field 
theory and are therefore natural physical quantities. In particular, if one of these poles 
comes close to the real frequency axis then we may be able to think of it as an emergent 
quasiparticle excitation of sorts (although it may not have a finite residue [H]). It may 
therefore be sensible to think about the isolated individual contribution of this pole to 
physical quantities. In our computation it was precisely such poles near the real axis which 
lead to the de Haas- van Alphen oscillations. More generally it might be interesting to 
compute the one loop contribution of such poles to other quantities such as the electrical 
conductivity. 

We primarily examined the low temperature regime in this paper, as the sum over 
quasinormal modes simplified in this limit and the nonanalytic motion of a specific pole 
contained the physics of interest. In general one would like to exhibit quantum oscillations 
at the higher temperatures more relevant for comparison with recent experiments. In par- 
ticular, it would be interesting to obtain a strong coupling analogue of the Lifshitz-Kosevich 
formula . To deal with higher temperatures it will likely be necessary to Poisson resum 
the Landau levels. Before this can be done, the asymptotic behavior of the quasinormal 
poles will need to be well characterised. We took some first steps in this direction above by 
computing the WKB form of the quasinormal poles at a fixed Landau level i. It would also 
be of interest to study the charged fermionic quasinormal modes numerically in a similar 
way to our studies of charged bosons in this paper. 

At the low temperatures we have studied, quantum oscillations are usually closely con- 
nected to the quantum Hall effect. It seems possible that the novel nonanalytic behavior 
of the low temperature free energy as a function of the magnetic field in theories with 
gravitational duals may indicate nonconventional quantum Hall states. 

In our numerical studies of the quasinormal frequencies of charged bosons we uncovered 
some nontrivial dependence of the pole locations on the charge of the field and the back- 
ground magnetic field. It would be very interesting to improve the numerical determination 
of the positions of the poles, perhaps by implementing a better numerical algorithm than 
the rather general method used in this paper. A more accurate knowledge of the 'pole 
dancing' will be necessary to establish unambiguously whether there are nonanalyticities 
as a function of the magnetic field due to bosonic quasinormal modes 'disappearing' into 
branch cuts. Also, it would be useful to be able to extend the numerics further into the 
lower half frequency plane and determine whether the branch cut terminates or whether it 
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continues asymptotically. 
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A Schrodinger form 



The equation ( 74 ) can be converted into a Schrodinger form (albeit with a complex potential 
at general z) 

-^ + V{r,)^ = L^X^, (102) 

by setting 

'^=717^^' = ^2- (103) 

This shows that the required norm is 

I dr,\^\^ = J ^|c/>|2<cx). (104) 
The potential is, written in terms of r, 

V = L m + Kj^r [ z — I 1 ~a ~a TTTT • (105) 



/V V r+yy 4 4 4 16/ 

At the boundary at infinity, r = 0, the potential tends to V^(0) = |-|-(Lm)^ which is positive 
so long as the scalar field satisfied the Breitenlohner-Freedman bound. This naturally 
suggests a continuum of eigenvalues at A > jp- + -m?. 

For self-adjoint operators one can prove that a complete basis of eigenstates lives in 
an extended Hilbert space which allows certain non-normalisable functions. At vanishing 
charge (q = 0) and pure imaginary frequencies z = icOn the operator M{z, £) is self-adjoint 
and positive. From the definition of the extended Hilbert space (see e.g. (5^) one can easily 
check that in this case there is a basis of eigenfunctions for which the condition A > +m? 
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holds. While the operator is not self-adjoint on the imaginary z axis with nonzero charge 
and chemical potential, these do not alter the boundary behavior of the field and we expect 
the continuum to persist as a basis of states. This statement is routinely used for charged 
fields in flat space and finite chemical potential. The only generalisation we are making is 
to consider a curved spacetime background. 



B Damped harmonic oscillator 



In this appendix we show how some of the methods we have developed in this paper can be 
applied to the very simple case of a damped harmonic oscillator. 
Consider the retarded Green's function 

1 



G\z) 



(106) 



-(z — /x)2 — 2i7z + m?- ' 
with < m and 7 > 0. This is just the response function of a damped simple harmonic 
oscillator. All the poles of G^{z) are in the lower half frequency plane. Associated with 
this retarded Green's function, we have the Matsubara Green's function defined on the 
imaginary frequency axis by 

1 



We now use the spectral representation to write 



G(z) 
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(107) 



(108) 



Note that G{z) = G^{z) in the upper half plane, but G{z) 7^ G^{z) in the lower half plane. 



However, on the imaginary frequency axis, with z = itOn, the expressions in Eq. (108) and 



(107) agree for all positive and negative u>n- Note that (108) defines an analytic function 



G{z) which has a branch-cut on the real z axis, and no poles anywhere. 

We can use the spectral representation to perform the frequency summation 



, . -c 



TT 



Im 



1 



1 



This object is of interest because it is related to the free energy by 

dJ^ 



dw? 



TY,G{iuJn). 



We can try to simplify (109) by using the position of the poles of G^: 

Ci 



(109) 



(110) 



:ill) 
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with Im[zj] < 0. Explicitly the poles are 2:1^2 and they obey 



Zl + Z2 
Z1Z2 
Zl,2 



2fi — 2i7 , 

2 2 
fi — m , 



fll2) 



Note that at large z the left hand side of (111) goes like ~ —Xjz^ which implies from the 
right hand side that 

^Ci = 0. (113) 

i 

This condition may often hold more generally independently of the number of poles, and 



was implicitly used in performing the frequency summation in Eq. (109). So we have 

r & 1 1 



.1 — rv-i '* 



Convergence of the integral at $7 = requires that 



Zi — $7 



114) 



J^Im 



0, 



(115) 



which is also seen to hold from Eqs. (Ill) and (106). Again, we expect Eq. (115) to hold 



more generally. It does not seem to be possible to simplify the integral in Eq. (114) further, 
in general. However, at T = 0, we obtain 



T ^ G{iujn) = - ^ Im \ci log(zj) 

Now we turn to the free energy. Let us write 

[G^(z)]-i = /(z)+m2, 



(116) 



Then, near a pole zi 



But f{zi) = —m? and so 



and hence 



So we have 



G^{z) 



f'iZi 



1 1 

f'{Zi) {Zi - z) 



dzi 



dm? 



dm? {zi — z) 



dzi 



Ci 



* dm^ 
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(117) 

(118) 
(119) 
(120) 
(121) 



This is an interesting result linking the residue and location of the 'quasinormal poles'. 



With this expression, we can write the constraints in Eq. (113) and (115) as 

d 
dm? 



0, 



d 

dm? 



Im l^log Zjj = . 



For the frequency summation, we have 

d^ d 



/OO JQ ^ 



1 



,n/T _ I 



Integrating this, using ( |110[ ), we obtain 

dn 



(122) 



(123) 



(124) 



up to m-independent terms. This expression is analogous to (56) in the main text. 
For the damped harmonic oscillator, we obtain as T — > 0: 

f°° dn^ „ , ^,,4 7 i"^ ndn 

^iT^o = — Im log(zi + --^^ / -Wt — ^ 

- - Im [zi log Zi] 



3 (m^ — iJ?) ' 



(125) 



where the Zi are given in Eq. (112) 



Alternatively, we can write 



T y -^^ + T y 



a;„>0 



Zl - ILOn 



(126) 



Integrating this with respect to m? we obtain 



T y^ \og{zi - iujr 



)+Ty\og{zi 



(127) 



This formula is analogous to (54) in the main text We can now evaluate the frequency 



summation by using the identities of zeta function regularization in equation (55). This 
yields 



^ = E 

i 

= E 



n=l 

OO 



n=0 

OO 



IZ^ 

2ttT 



r E >°4» + S) + ^ E - 1^) - 

. n=0 ^ ^ n=0 ^ ^ 



2ttT 



-ryiog 



\zi\\T{iz^/i2TTT))\' 



(128) 
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This result is analogous to (52) in the main text. 

As T ^ at fixed Zi, we can use the asymptotic expansion of the T function 



logr(x) 



(129) 



Working with all the terms in Eq. (129), we have the complete expression as T ^ 0, 



assuming all the Zi remain finite in this limit 



1, 



-Im 



vr 



1 T 
Zilog{izi)^ + - log 



27rr 



Zi 



log(l/(27rer)) 



^Im(z,) + ^ J]Re(zi) 

i i 

^ ^ Im [z, log z,]-^Y. I™ (^) + ^(^')- 



(130) 



The above expression is very general. Specializing to the values in Eq. (112) we obtain 

2tt jT"^ 



3 (m^ — fi"^) 



+ 0(r^)(131) 



Apart from the first two m-independent terms, this matches precisely with Eq. (125). 



C Asymptotic frequencies 



In this appendix we outline the derivation of the WKB formula (100) for the charge quasi 



normal modes in the main text. We will follow the notation of (STl which uses different 



coordinates to those in this paper. Let us therefore rewrite our equation (75) as 

2 



d 



F d 



dR \ R'^ dR 



F 



where 



F{R) = 1-4 



UJ 



R 



qfi 1 



R 



+ (KeR^ + 



<^>{R) = 0, (132) 



^1 - 7rr+T) + 



(3 - 47rr+T) . 



(133) 



R 

The boundary is at ii = and the horizon is at = r_|_ . We are interested in the quasinormal 
frequencies of this equation. 



We want to transform equation (132) into the notation of Ref. So we define 



f{r) 
4>{r) 



R ' 

r+R-^F{R) 



r 
r+ 

r^{R) 



4(1 - vrr+T) (3 - Airr+T) 



+ 



(134) 
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Now the boundary is at r = cxd and the horizon is at r = 1, and the differential equation 
transforms to 

d [ M\ [m^ Ki,r+ 1 df\ ( 



-f-r(f^] + ( — + ^ + -f]f<P=(^-<ii^- -] <P- (135) 

ar \ ar J \rj^ r ar J \ ^ J 

This is similar to Eq. (3.2) of Ref. |51], allowing us to map some of their results to ours. 
We will need the roots of /(r) below. The roots are at r = with 

n = 1 ; /'(ri)=47rT, 
r2 = r_/r+ ; < r2 < 1 ; /'(r2)<0, 
rsA = -(l+r_/r+)/2±iCT. (136) 

where a is real. These roots satisfy some useful identities: 

4 



p=i 

tfh-- ''''' 

The WKB formulae will involve contour integrals which are most easily expressed in 
terms of a coordinate x, related to r by 

dv 

dx = Tr^- (138) 
fir) 



This can be integrated to obtain 



4 

^(0 = >^777^1og( 1- - ) • (139) 



The WKB limit of the quasinormal modes is captured by matching together the behavior 



of solutions to ( 135 ) near several special values of r. We now consider these values in turn. 



(black hole singularity) 



The branch-cuts in Eq. ( |139 ) should be chosen so that there is no monodromy around r = 0, 



and also none around r = oo. Then the expansion near r = is 



+ . . . , (140) 



3(3 - 47rr+r) 
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and the differential equation (135) becomes 



-(j) = UJ (p. 



This has the same form as that in Section 3.3.2 of Ref. [51] with 

j = 1/3. 



fl411 



(142) 



Note that for |x| <C 1, the term proportional to q in Eq. (135) is subdominant to either 
the or the l/x"^ terms in Eq. (141 ). This means that we can solve Eq. (141 ) in terms of 



Bessel functions. In the regime where ojx ^ 1, but with |x| <C 1, we have the solution as in 
Eq. (3.40) of |5l]: 



(x) ~ 2i?_|_ cos {iox — a+) + 2i?_ cos {ivx — a_) 



(143) 



where = 7r(l it j)/4. 

However, for the matching below, we have to extend this solution to x ~ 1. In this 



regime we can use the WKB method, where the right hand side of Eq. ( 135 ) dominates over 



the potential on the left hand side. From this method, with ^ 1, we find that we must 

(r-1) 



replace Eq. ( 143 ) by 



(x) ~ 2Ba. cos 



dx I to — qjjL- 
V r 



a+ +2B^ cos 



dx [ uj — qfi- 
V r 



a- 



(144) 



The integral is to be taken over a suitable contour in the complex plane, which we will 
discuss below. 

r = oo (asymptotic AdS region) 

Here we have 



Xq 



x{r — > oo) 
1 



4 



log 



(145) 



where xq is a complex number which will play an important role below. Despite the am- 
biguity in the branch cuts of the logarithms, the value of xq is unique: the identities in 



Eq. (137) ensure that there is no monodromy around r = oo provided none of the branch 



cuts extend to r = oo. For x close to xq, we have 



Xq — X 
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(146) 



and the differential equation Eq. ( 135 ) becomes 

w? + 2 



+ 



(p = LO (p . 



(147) 



dx^ (xo — x)2 

Again, the qfj, term is subdominant everywhere for |x — xq] <C 1. In the limit of large uj this 
is the same as in Section 3.3.2 of Ref. [5T] with 



Jc 



\/4m2 + 9. 



As X — > Xo, the solutions of Eq. (147) are of the form 

,(-l±V4m2+9)/2 



(148) 



(149) 



implying that in terms of our original variables ^> ~ 

^(3±v^W+9)/2 

as expected. 

The connection to the regime where |x — xo| <C 1, but uj{xq — x) ^ 1, can be performed 
using Bessel functions, as in Ref. |51j. Requiring falloff at infinity yields the result above 
(3.41) in m- 

(/>(x) ~ 2C+ cos(a;(x - Xo) + /3+) , (150) 



with /3+ = 7r(l +ioo)/4. As for Eq. (143), we have to extend this result for xq — x ~ 1, but 



\uj\ ^ 1. Here the WKB method shows that Eq. (150) is replaced by 



(x) ~ 2C+ cos ( / dx [ UJ — qfi- ^ ) + j , 



where the contour of integration is left unspecified for now. 



r = 1 (black hole horizon) 



:i5ii 



Near the horizon of the black hole 



The equation (135) becomes 



and so the solution is 



log(r — 1) 
X — — > — oo , 



(r) « (r - l)±W(4^r) ^ 



(152) 



(153) 



(154) 



which is as expected. We will want to impose ingoing boundary conditions, as in the main 
text. 
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Asymptotic frequencies 

We now see that the determination of the frequencies of the quasi-normal modes can be 
mapped onto the solution in Section 3.3.2 of Ref. [Slj, after using the values of j, joo and 
xq quoted above, and including the phase shifts from the qfi term. 



Matching the behavior near r = in Eq. (144), with that near r = oo in Eq. (151 ), we 



find that the condition (3.41) in [SIJ is replaced by 

exp^2.(/ dx{u^-qf,^—^]-P+]\ = ^^±—^^—^, (155) 







where the integral is to be taken along contour B connecting these points in Fig. 15 of 
Ref. [51] • Namely, it is the Stokes line satisfying Iniux = 0. This is crucial in order to be 
able to accurately distinguish ingoing and outgoing modes at the horizon. 

The matching between r = and the horizon r = 1 remains as in Ref. |51j . and we have 
the result above their (3.42). 

^^g-3ia+ _^ 5_e-3*"- = 0. (156) 

The final result for uj is obtained by solving these last two equations, and this shows 
that Eq. (3.42) of Ref. [51j is replaced by 

^0 / (r-l)\ 1, / /vr\\ 7r(2 + V4m2 + 9) , , 

dxiio-qfi^ '- =n7r+-log(2cos(-)) + ^ ^ ^ (157) 







2f°V V6// 4 



where n is an integer, and the integral is along contour B. This is equation (100) in the 
main text. In the main text we transformed back to our coordinates. 
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